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The theory of a flux steady-state (avalanche) formation is presented for the simplest model of a 
real sand pile within the framework of Lorenz approach. The stationary values of sand velocity and 
sand pile slope are derived as functions of controlling parameter (externally driven sandpile slope) . 
The additive noises of above values are taken into account to build the phase diagram, where the 
noise intensities determine a domain of the avalanche appearance. This domain shows to be crucial 
to the noise intensity of vertical component of sand velocity. 
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I. INTRODUCTION 

In recent years considerable study has been given to the theory of self-organized criticality (SOC) [1] that ex- 

3 . plains spontaneous (avalanche-type) dynamics, unlike the typical phase transitions that occur only when a controlling 

parameter is tuned to a critical value. A main feature of the systems displaying SOC is that they are distributed 

over avalanche sizes, so that SOC models are mostly studied by making use of scaling-type arguments supplemented 

with extensive computer simulations (see [2]). On the contrary, in this paper we put forward the theory of a single 

Q avalanche formation. 

The SOC behavior appears in a variety of systems such as a real sand pile (ensemble of grains of sand moving 
along increasingly tilted surface), intermittency in biological evolution [3], earthquakes and propagation of forest- 
fires, depinning transitions in random medium and so on [4]. Among the above models the simplest and most widely 
studied, analytically [5], [6] and numerically [1], [7], [8], are the sandpile models. The field theory [9], based on 
nonlinear diffusion equation, has failed to account for the basic feature of self-organized systems - their avalanche 
dynamics. The reason is that there is no feedback between open subsystem and thermostat within the framework of 
the one-parameter approach employed in [9]. Recently the two-parameter theories were set forth in [10,11], where 
Q\ • the thermostat degree of freedom was either controlling parameter [10], or conjugate field [11]. In the mean field 
approximation, the approach of [11] allows to obtain critical exponents governing scaling behavior of self-organized 
system. Our approach is to take into consideration the complete set of degrees of freedom. Owing to this, not only the 
complimentary results of [10,11] are reproduced, but it is possible to obtain the self-consistent analytical description 
of the single avalanche formation. 

This paper is organized as follows. In Sec. 2 the self-consistent theory of the formation of a flux steady-state is 
presented. It enables us to treat the problem on the basis of the unified analytical approach. Sec. 3 deals with 
accounting additive noises of the sand velocity components and sandpile slope. It is shown that the noise intensities 
increase the possibility of emergence of the SOC regime. Sec. 4 contains the obtained results. 



(N 



O 



Ch . II. NOISELESS CASE 

o 



X 



Within framework of the simplest model of a sandpile, a dependence y — y(t, x) defines its surface at given instant of 
time t. Locally the flow of sand can be described in terms of three quantities: the horizontal and vertical components 
of the sand velocity, x = dx/dt, y — dy/dt, and the surface slope y' — dy/dx. The key point of our approach is 
that the above degrees of freedom are assumed to be of dissipative type, so that, when they are not coupled, their 
5« \ relaxation to the steady state is governed by the Debye-type equations: 

dx x , . 
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where t x , r y and t$ are the relaxation times of the velocity components and the slope, respectively. Eqs. dTa) - 
(|lq) imply the sand is at rest in the stationary state, x — y = and the equilibrium slope y' = y' a ^ plays the role 
of a controlling parameter. 

Since the motion of sand grain along different directions is not independent, Eq.(|la|) should be changed by adding 
the term / = y/j due to liquid friction force along the y-axis (7 is the kinetic coefficient). Then, we have 

t x x — — x + a~ l y (2) 

where a = j/t x . Note that, owing to the diffusion equation y — Dy" (D is the diffusion coefficient), the friction force 
appears to be driven by the curvature of the sandpile surface: 

/ = {Dh)v". (3) 

On the other hand, when x = (stationary state) , solution of Eq. (0) defines the tangent line y — ax + const, so that 
the friction force / — t x 1 x is proportional to the horizontal component of the sand velocity. Taking into consideration 
the relation (||) and using the chain rule dy' /dt = if + y"x, from Eq.(fTq) one obtains the equation of motion for the 
slope: 

T S y' = (y' -y')-(Ts/D)yx. (4) 

Following the same line, the equation for the vertical component of the velocity can be deduced 

'.»--* + **■ T v ^(^~)- < 5 > 

Note the higher order terms are disregarded in Eq.(g) and the renormalized relaxation time T y depending on the 
stationary slope Wq is introduced. 

Eqs.(||), (||), @ constitute the basis for self-consistent description of the sand flow on the surface with the slope 
y' driven by the control parameter y' . The distinguishing feature of these equations is that nonlinear terms that 
enter Eqs.Q, (O) are of opposite signs, while Eq.(g) is linear. Physically, the latter means that on the early stage 
the avalanche begins moving along the tangent y = ax + const. The negative sign of the last term in Eq.(0) can be 
regarded as a manifestation of Le Chatelier principle, i.e. since an increase in the slope results in the formation of an 
avalanche, the velocity components x and y tend to impede the growth of the slope. The positive feedback between 
x and y' in Eq.rtq) plays an important part in the problem. As we shall see later, it is precisely the reason behind the 
self-organization that brings about the avalanche generation. 

After suitable rescaling, Eqs.(0), (0), (H) can be rewritten in the form of the well-known Lorenz system: 

u = —u + v, (6a) 

ev = —v + uS, (6b) 

6 S=(S -S)- uv, (6c) 

where u = (Ty/r^irs/D) 1 / 2 ^ v = {r v /T x ) l l 2 (r s ID) 1 /' 1 y/a^udS 

{Ty/Tx) 1 ' 2 ^ /a are the dimensionless velocity components and the slope, respectively; e = t v /t x , S = tsJt x and 
the dot now stands for the derivatives with respect to the dimensionless time t/r x . In general, the system (|a]) - ( pq ) 
cannot be solved analytically, but in the simplest case, where e <C land 5 ^C 1, the vertical velocity v and the slope 
S can be eliminated by making use of the adiabatic approximation that implies neglecting of the left hand sides of 
Eqs.(pq), (pq). As a result, the dependencies of S and v on the horizontal velocity u are given by 

a S S u 

S = -, ^j v = "; 7- (7) 

1 + u 2 l + u 2 y ' 

Note that, under u is in the physically meaningful range between and 1, the slope is a monotonically decreasing 
function of u, whereas the velocity v increases with u (at u > 1 we have dv/du < and the flow of the sand becomes 
turbulent). 

Substituting second Eq.(Q) into Eq.(|6a|) yields the Landau-Khalatnikov equation: 



dE 
du 



(8) 



with the kinetic energy given by 



E=~u 2 -±S ln(l + u 2 ). (9) 



For So < 1, the it-dependence of E is monotonically increasing and the only stationary value of u equals zero, uq = 0, 
so that there is no avalanches in this case. If the slope So exceeds the critical value, S c = 1, the kinetic energy assumes 
the minimum with non-zero steady state velocity components u e = v e = (So — l) 1 ' 2 and the slope S e = 1. 

The above scenario represents supercritical regime of the avalanche formation and corresponds to the second-order 
phase transition. The latter can be easily seen from the expansion of the kinetic energy (Q) in power series of u 2 <C 1. 
So the critical exponents are identical to those obtained within the framework of the mean field theory [11]. 

The drawback of the outlined approach is that it fails to account for the subcritical regime of the self-organization 
that is the reason for the appearance of avalanches and analogous to the first-order phase transition, rather than the 
second-order one. So one has to modify the above theory by taking the assumption that the effective relaxation time 
t x (x) increases with velocity u from value r x (l +m) , m > to t x [12]. The simplest two-parameter approximation 
is 

m 

1 + t-t-, r, ( 10 ) 



t x (u) l + {u/u ) 2 

where < Uo < 1. The expression for the kinetic energy (9h then changes by adding the term 



m _ 2 . 




AE=-utln[l+[-) ) (11) 



and the stationary values of u are 



u? = uoo (l T [1 + uluooHSo - S c )] 1/2 ) V2 , (12) 

2«q = (So — 1) + S c ul, S c = 1 + m. 



The upper sign in the right hand side of Eq. ( p. 2D is for the value at the unstable state u m where the kinetic energy 
E + AE has the maximum, the lower one corresponds to the stable state u e . The corresponding value of the stationary 
slope 



1 + u 2 00 + \/ (1 + u 2 0Q f - (1 - u 2 ) So 



s m = — v 1 : 2 — (is) 



i-«S 



smoothly increases from the value 



= l + uoyJm/(l-v$) (14) 



at the parameter So = S c o with 

ScO = (l - U§) S'min (15) 

to the marginal value S c — 1 + m at So — S c . The ^o-dependencies of u e , u m , andS'e are presented in Fig.l. As is 
shown in Fig. la, under the adiabatic condition ts ^C t x is met and the parameter ^o slowly increases being below 
S c (So < S c ), no avalanches can form. At the point So = S c the velocity u e jumps upward to the value \/2uoo 
and its further smooth increase is determined by Eq.(n2l). If the parameter So then goes downward the velocity u e 
continuously decreases up to the point, where So = S c q andu e = uoo- At this point the velocity jump-like goes down 
to zero. Referring to Fig. lb, the stationary slope S e shows a linear increase from to S c with the parameter So being 
in the same interval and, after the jump down to the value (1 — it 2 ,) -1 a l ^o = S c , S e smoothly decays to 1 at So 3> S c . 
Under the parameter So then decreases from above S c down to S c o the slope grows. When the point Jig ) is reached, 
the avalanche stops, so that the slope undergoes the jump from the value ([14]) up to the one defined by Eq. (fl5|) . For 
So < Sco again the parameter Se does not differ from So- Note that this subcritical regime is realized provided the 
parameter m, that enters the dispersion law dlQ), is greater than 



m m in = 2 - (16) 

1 U Q 

Clearly, according to the picture described, the avalanche generation is characterized by the well pronounced 
hysteresis, when the grains of sand initially being at rest begin to move downhill only if the slope of the surface 
exceeded its limiting value S c = 1 + m, whereas the slope S c o needed to stop the avalanche is less than S c (see 
Eqs.(O), (|l5|)). This is the case in the limit ts/t x — ► and the hysteresis loop shrinks with the growth of the 
adiabaticity parameter S = ts/t x . In addition to the smallness of 5, the adiabatic approximation implies the ratio 
t v /t x = e is also small. In contrast to the former, the latter does not seem to be realistic for the system under 
consideration, where, in general, r y rj t x . So it is of interest to study to what extent the finite value of e could change 
the results. 

Owing to the condition S <C 1, Eq.(J6q) is still algebraic and S can be expressed in terms of u and u. As a result, we 
derive the system of two nonlinear differential equations that can be studied by the phase portrait method [12]. The 
phase portraits for various values of e are displayed in Fig. 2, where the point O represents the stationary state and 
the point S is related to the maximum of the kinetic energy. As is seen from the figure, independently of e, there is the 
universal section that attracts all phase trajectories and its structure is appeared to be almost insensitive to changes 
in e. Analysis of time dependencies v(t) and u(t) reveals that the velocity components slow down appreciably on this 
section in comparison to the rest parts of trajectories that are almost rectilinear (it is not difficult to see that this 
effect is caused by the smallness of 5) . Since the most of time the system is in vicinity of this universal section, we 
arrive at the conclusion that finite values of e do not affect qualitatively the above results obtained in the adiabatic 
approximation. 

III. NOISE INFLUENCE 

To take into account additive noises of the velocity components u, v and the slope S it needs to add to right-hand 

-V 
S 



parts of Eqs.(6a) - (pel) the stochastic terms I u £, ij £, ij £, respectively (here the noise intensities I u , v ,s are 



measured in units {t x ~Jt v ){D/ts), a 2 {t x / t v )(D / ts) 1 a 2 (T x /r y ), correspondingly, and £(£) is (^-correlated stochastic 
function) [13]. Then, within the adiabatic approximation Eq.(0) acquires the stochastic addition 



{ll /2 +[ll /2 g v {u) + iy 2 g s {u)\)m 



where we introduce the multiplicative functions gs(u) = ug v (u) = u/(l + u 2 ). As a result the extreme points of the 
distribution P{u) oc exp{— U(u)} of the stochastic variable u is given by the effective energy [14] 

U(u)=lnl(u)+ f^/^du (17) 

J I(u) 

where the bare energy E is determined by Eq.(H) and the expression for the effective noise intensity 

I{u) = I u + I v g 2 ,{u) + I s g 2 s (u) (18) 

follows from the known property of additivity of squares of variances of independent Gaussian random quantities 
[13]. Combining expressions (g), (|l7|), J18]), we can find the explicit form of U(u), which is too cumbersome to be 
reproduced here. Much simpler is the equation 

x 3 -S Q x 2 -2I s x + 4:(I s -I v ) =0, x=l + u 2 , (19) 

which defines the locations of the maxima of distribution function P(u). According to Eq.(|19|), they are insensitive 
to changes in the intensity of noise I u of the velocity component u, but are determined by the stationary value So of 
the sandpile slope and the intensities I v , Is of the noises of vertical velocity component v and slope S, which acquire 
the multiplicative character in Eq.fllgj). Hence, it can put for simplicity /„ = and Eqs.(|9|), (|L7j), ( |l8| ) give the follow 
expression for the effective energy: 

U(u 



1 

2 

(1 -i)(l-S - i) Hi + u 2 )} + I S InfoKu) + igl{u% i = I v /I s . 



- + (2 - S - i)u 2 + (20) 



According to Eq.(|l^), the effective energy (|2^) has a minimum at u — if the stationary slope So does not exceed 
the critical level 



S c = l + 2/ s ~4/„, 



(21) 



whose value increases at increasing intensity of noise of the sandpile slope, but decreases with one of the velocity. In 
this case, sand grains are at rest. In the simple case I v = 0, the avalanche creation corresponds to solutions 



1 



(22) 



So - 3 + v /(3-S ) 2 +4(2S'o-3 + 2/s) 
which are obtained from Eq.(n9) after eliminating the root u 2 — 0. The magnitude of this solution has its minimum 

< = \ [(So - 3) - y/(S + 7)(S -lj\ (23) 

on the line defined by expression ( UW with /„ =0. At So < 4/3, the roots ±u c are complex, at So — 4/3 they become 
zero, at So > 4/3 they are real, and u + = — u_. In this way, the tricritical point 



So = 4/3, Is = 1/6 



(24) 



corresponds to the appearance of roots u± ^ of Eq.([l^) (avalanche creation). If condition ( j2l| ) is satisfied, the root 
u = corresponds to the minimum of the effective energy ( |20|) at So < 4/3, whereas at So > 4/3 this root corresponds 
to the maximum, and the roots u± to symmetrical minima. 

Now we find another condition of stability of roots u±. Setting the discriminant of Eq.(|19|) equal to zero, we get 
the equations 



Is = 0, 



Is -Is 



27 



-r 1-^ -^ 



st 



1 = 0, 



the second of which gives 



21* = 



^!fi_^V^ 



± 



!Z(i_^V^ 2 



1/2 



-2S n J 



(25) 



(26) 



This equation defines a bell-shaped curve So(Is), which intersects with the horizontal axis at the points Is 
Is = 27/2, and has a maximum at 



So — 2, 



la = 2. 



and 



(27) 



It is easy to see that for /.„ = this line touches the curve (|2l]) at point (|24[). 

Let us now consider the more general case of two multiplicative noises I v ,Is 7^ 0. Introducing the parameter 
a = 1 — i, i = I v /Is and the renormalized variables / = Is/a 2 , So = Sq /a, u 2 = (1 + u 2 )/a — 1, at i < 1 we may 
reproduce all above expressions with the generalized energy U/I in Eq.ptl). Then the action of the noise of the 
vertical component of the velocity v is reduced to the renormalization of the extremum value of the horizontal one by 
the quantity (a -1 — l) 1 / 2 , so that the region of divergence u « becomes inaccessible. 

The condition of extremum of the generalized energy (20) splits into two equations, one of which is simply u = 0, 
and the other is given by Eq.([19|). As pointed out above, analysis of the latter indicates that the line of existence of 
the zero solution is defined by an expression (Elf). The tricritical point has the coordinates 



So 



;(1 



(1 + 81,,). 



(28) 



The phase diagram for the fixed intensities /„ is shown in Fig. 3. Here the curves 1, 2 define the thresholds of absolute 
loss of stability for the fluxless and flux steady-states, respectively. Above line 1 the system occurs in a stable flux 
state, below curve 2 it is in fluxless one, and between these lines the two-phase domain is realized. For I v < 1/4 
situation is generally the same as in the simple case I v = (see Fig. 3a). At /„ > 1/4 the SOC regime is possible for 
small intensities Is of the slope noise (Fig. 3b). According to ( p8| ) the tricritical point occurs on the Is axis at I v = 1, 
and for the noise intensity /„ larger than the critical value I v = 2 the stable fluxless state disappears (see Fig. 3c). 



IV. CONCLUSION 

According to the above consideration, the dissipative dynamic of grains flow in a real sand pile can be represented 
within the framework of Lorenz model, where the horizontal and vertical velocity components play a role of an 
order parameter and its conjugate field, respectively, and the sandpile slope is a controlling parameter. In Sec. 2, the 
noiseless case is examined to show that an avalanche creates if the externally driven sandpile slope y' is larger than 
the critical magnitude 7(T a; Tj,)~ 1 / 2 . In this sense, the systems with small values of the kinetic coefficient 7 and large 
relaxation times t Xi t v of the velocity components are preferred. However, the sand flow appears here as a phase 
transition because the spontaneous avalanche creation is impossible in the noiseless case. Taking into account the 
additive noises of the above degrees of freedom, we show in Sec. 3 that the stochasticity influence is non-essential for 
the horizontal velocity component and is crucial for the vertical one. The SOC appears if the noise intensity of the 
latter exceeds the value 2^ 2 (D-y 2 /t x t v ts), provided the noise is small for the sandpile slope. If the noise intensity of 
the vertical velocity component is larger than the critical value 2(Dj 2 /t x t v ts), the fluxless steady-state disappears 
at all. 
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FIGURE CAPTIONS 

Fig.l. The So-dependencies of a) the velocities u e , u m , and b) the equilibrium slope S e . The arrows indicate the 
hysteresis loop. 

Fig. 2. Phase portraits in the v — u plane at So = 1.25S* c o for a) e = 1CP 2 ; b) e = 1; c) e = 10 2 . 

Fig. 3. Phase diagrams for fixed values /„ of the noise intensities of the vertical velocity component: a) /„ = 0, 
b) I v = 1, c) I v = 2. Curves 1 and 2 define the boundary of stability of avalanche and non-avalanche phases; A - 
avalanche phase, N - non-avalanche. 
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